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The masses of the Ds(O^) mesons are investigated from a view-point of ordinary light-heavy 
system in the framework of the Gaussian sum rules, which are worked out by means of the Laplacian 
transformation to the usual Borel sum rules. Using the standard input of QCD non-perturbative 
parameters, the corresponding mass spectra and couplings of the currents to the Ds{0 ) mesons are 



obtained. Our results are m 



D,(0-) 



1.968±0.016±0.003 GeV and m 



DAO+) 



2.320±0.014±0.003 



GeV, which are in accordance well with the experimental data, 1.969 GeV and 2.317 GeV. 
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Recently, the experimental observations of D meson 
and the corresponding theoretical manipulations attract 
much attention in the research of particle physics [l|-llO|. 
Many Quantum Chromodynamics (QCD) calculations 
about D mesons are in good agreement with the exper- 
imental data within the theoretical uncertainty except 
for the unexpected low mass of £'o+(2317) |lll - [l4 |. In 
Ref. [H, the masses of the excited (0+, 1+) and (1+, 2+) 
doublets for the cs system are calculated to the l/rric or- 
der in the sum rules based on heavy quark effective the- 
ory (HQET) , and the D^(2317) and Ds{24:60) are iden- 
tified as the (0+, 1+) doublet so that the mass splitting 
in this doublet is well reproduced. We note that most 
of these theoretical calculations are based on QCD Borel 
sum rule (BSR) which emphasize the contributions of 
lowest resonance state, and have shown the power in 
the investigation of the non-perturbative properties of 
hadron bound states. However, excited states and con- 
tinuous spectrum can produce background interference. 
Generally, this background interference increases with 
the mass of the hadron state considered, and this restrict 
the application range of BSRs. It is noticed that the 
QCD Gaussian sum rule (GSR) developed later empha- 
sizes only the contribution of the hadron state considered 
as seen from the appearance of the Gaussian distribution 
function, and has more clean background in comparison 
with the BRS. From this reason, the GSR may, in princi- 
ple, work better than the BSR. At least, both sum rules 
should give almost the same results because they are de- 
rived from the same underlying dynamical theory. 

On the other hand, whether the Dg meson's spec- 
troscopy observed by experiments can be derived out 
from the elementary theory of the strong interactions 
is a key point for testing the QCD non-perturbative 
mechanism, and serves also to be the starting point of 
calculating the important physical processes, such as 
the Ds meson decays and etc, from the first princi- 
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pie. However, the results obtained by using the BSRs, 



GeV 

data. 



4^= 2.48 ± 0.03 GeV and m£,^(n-) = 1.94 ± 0.03 



li|, are obvious inconsistent witji the experiment 



^040+) 



2.317 ±0.0013 Ge 



and m 



Ds{o- 



1.969 GeV. For checking the correctness of QCD and ex- 
ploring the non-perturbative mechanism, it's necessary 
to do a further study. 

In order to make a cross check between various types 
of QCD sum rules, namely the BSRs[ll| and the HQET 
sum rules Is], we calculate the masses rriD^ of Dg mesons 
and the couplings fu^ of the corresponding currents to 
these mesons by using GSRs. 

Consider the gauge-invariant and Lorentz-covariant 
two-point function of the current Ji (x) corresponding to a 
resonance with the quantum numbers 1 = 0^,0"*" |16l . |l7| 



n.(g^) 



d^ 



3'^-(0|rj.(x)4(0)|0), (1) 



where |0) is nonperturbative QCD vacuum, and Ji{x) are 
the pseudoscalar and scalar currents 



Jo-(x) = iq(x)jr,c{x), 
Jo+{x) = q{x)c{x), 



(2) 
(3) 



with q{x) and c(x) being the strange-quark and charm- 
quark fields at the point x, respectively. For the invariant 
functions H(g^) = Hi(g^)(i = P, S), we have a dispersion 
relations without any subtraction 



n(g2 



ImH(s) 

> '■ 

s — q-^ -\- % 



(4) 



Via this dispersion relation, the QCD sum rule (QSR) 
was constructed by equating the contribution of opera- 
tor product expansion (OPE) and the phenomenological 
(PH) one related to the spectral function. The former is 
described as the product of Wilson coefficients and non- 
perturbative QCD vacuum condensates or quark masses, 
while the latter is parameterized by hadronic quantities 
such as resonance masses, couplings and the continuum 
threshold, etc. Thus, the QSR can be represented in a 
simple form. 



' 1 

dsW{s)- [imH 

TT 



PH/ 



s) - ImHOPE(s)] = 0, (5) 



where W{s) is an arbitrary weight function being ana- 
lytic except for the positive real axis starting from the 
lower mass squared, m^. According to the most success- 
ful application of sum rules to mesons and baryons 16], 
the phenomenological spectral function is assumed to be 
saturated by one narrow-width resonance and a contin- 
uum 

-ImnP«(s) =M(s~m|j) + -ImnOPE(s)0(s-So), (6) 

TT TT 

where sq is the QCD continuum threshold, the pole 
residue is of the form F — f]}rnj^ with /^ being the 
couplings of the lowest resonances with respective pari- 
ties to the hadronic currents and m^ a pole mass. The 
power k of mj^ in the pole residue is taken to match the 
maximum power of s in the asymptotic s-behavior of the 
spectral function. For s > sq, the hadronic continuum re- 
duces to the same form with that obtained by an analytic 
continuation of the OPE[il[, i.e. the perturbative terms, 
based on a hypothesis of the quark-hadron duality. 

Implementing the Borel transformation to the QSRs 
and performing the OPE at dimension d < 6 operators, 
the relations of the BSRs of the lowest 0^ cs-mesons are 
obtained to be Til 
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ms{ss)o 



8 \ TT /O 



(7) 



for the scalar 0+ and pseudoscalar 0^ channels respec- 
tively, where g denotes strong coupling constant, G^ — 
G/ij/G'"^ with Gf_ii, being the gluon field and a ■ G — 
Cjj.vG'^'^ , and {0)q represents the vacuum expectation 
value of a local composite operator 0(0) at the origin. 
Here, ?7is-correction is maintained at the first order [lll |. 
and the a^ correction to the perturbative contribution is 
given as the functions, Ro{x), as follows jlSl [l9| 



Roix) 



9 



+ 2Li2{x) + In a; ln(l - x) 



ln(l — x) + xln- 



2 X 

■ lux. 



(8) 



with the Spence function Li2{x) = — /q dtt~^ln{l — t). 
The running coupling constant as{s) appearing in the 
perturbative terms of Eq.(l7]) is approximated by a 
one-loop form [ll|, as{s) — 47r/[91n(s/Ag^^)] with 
^QCD — (0-25 GeV)^, which is determined to reproduce 
as{l GeV) ~ 0.5 ^. We note here that the contri- 
butions above the threshold sq on both sides of Eq.® 
are assumed to be equal to each other due to the quark- 
hadron duality at high scale, and have been removed out. 
In fact, our expression of Eq.([8]) is different from Ref.[ll], 
namely, some condensate contributions which remain fi- 
nite above sq is thrown off besides the perturbative con- 
tribution (s. the last line in Eq.([8|)). 

In order to derive the GSRs, we make the Laplacian 
transform on Eq.([7]) with the formula 21 
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-(s+s)(j 
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dse 



■-Imn(s) = G(-s,T). (9) 



This gives a procedure to construct G{—s, r) from Eqd?]), 
and hence G{s,t) by analytic continuation. The resul- 
tant expressions of GRSs for the Ds{0^) meson currents 
are 
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in the scalar 0^ and pseudoscalar 0^ channels, respec- 
tively. 

Taking the derivative with respect to s for both sides 



of Eq.dlOl), we get 

r2 2/2 -\ - 



.^fl. ,n, 



From Egs. dTUl) and pTj) . we obtain the sum rules for the 
masses of lowest resonances and couplings to the corre- 
sponding currents 



to± 



J.2 / - \ G{s,t) 
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(12) 
(13) 



In order to extract the values of physical quantities 
from the GSRs, we use the following standard values for 
QCD parameters appearing in the OPE^ll,J(see Tab|l|. 
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FIG. 1: The curves oiniR = rnjj^iQ-^ vs. r from GSRs, where 
s = 1.96^ GeV^ so^(o-) = 4.5 ~ 4.7 GeV^ 



TABLE I: QCD input parameters used in the analysis. 



Parameters 



References 



ms =0.11 ±0.01 GeV [16] 

rUc = 1.46 GeV [13] 

(nn>o = (-0.225 ± 0.025 GeV)^ [16| 

(ss)o = (0.8 ± 0.1) X (nn)o GeV^ [16| 

(^G^>o = (0.33 GeV)^ [18| 

(npo- ■ Gn)o = Afo (nn)o [22] 

Mo = 2 X (0.4 ± 0.1) GeV^ [22] 
{sga ■ Gs)o = Mo{ss)o 
as{nn)l = 0.162 x 10"^ GeV*' 
as{ss)l = (0.8 ±0.1)^ X as{nn)l GeV^ 
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For numerical calculation, we must determine the val- 
ues of s and the thresholds sq. To investigate the prop- 
erties of the considered resonance, the the value of s 
should approximately be set to be the corresponding 
mass squared, Tn|j , of the resonance. To suppress the 
continuum contribution, we require s < 7n|, . The condi- 
tions for determine the value of sq are: first, it should be 
grater than ?7i|j ; second, it should guarantee that there 
exists a sum rule window for our sum rules. We note that 
the upper limit of the sum rule window is determined by 
requiring the contribution of continuum to be lower than 
30% of the total, while the lower limit of that window 
is obtained by requiring the non-perturbative contribu- 
tions, proportional to a positive powers of cr, to be less 
that 30% (in fact, less than 10%) of the perturbative one. 
Therefore, the value of sq: the upper and lower limits of 
the corresponding sum rule window are determined in 
self-consistent manipulation. 

Then, using the GSRs of the mass and coupling con- 
stant of 1)5(0^), we can get the figures and numerical 
results shown below. Figs. [1] and [2] display the dependen- 
cies of the calculated masses m£)iQ±\ on r from GSRs, 
Figs. 13] and m plot the coupling fD,{o±) vs r. In Figs. 
[T] and m the thick horizontal-lines on the curves are the 



FIG. 2: The Gauss curves of tur = m,£,_,(o+) ^s. r, where 
s = 2.28^ GeV^, So,(o+) = 6.5 ■- 6.7 GeV^. 
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FIG. 3: The Gauss curves of Jb, = fDjo-\ vs. r, where 



s = 1.96^ GeV^ s 
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TABLE II: The numerical results of QCD sum rules. For 
comparison, we attach experimental avarage values observed 
[j-[^ and the masses of the first radial excitations predicted 
in Ref.[ll]. (GeV) 
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(GSR) 
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(exp.) 


(model) [12] 


0- 

0+ 


1.968 ±0.016 ±0.003 
2.320 ±0.014 ±0.003 


1.94 ±0.03 
2.48 ±0.03 


1.969 
2.317 
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3.067 
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0.158 ±0.006 ±0.003 
0.225 ±0.005 ±0.003 
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FIG. 4: The Gauss curves of /jj = 10^10+) ^^- '''i where 
s = 2.28^ GeV^ Sb^(o+) = 6.5 ~ 6.7 GeV^ 



Gauss stabilities of the sum rule windows. The length of 
the lines corresponds to the sizes of stable regions. 

In Ds{0~) channel with s = 1.96^ GeV , we cannot 
find any Gaussian stability below Sq- — 4.5 GeV^ and 
above sq- = 4.7 GeV , while between these two values 
the curves are stable. From these plateau regions, we get 
the resonance mass, mD^{o-), to be 1.968 ±0.016 ±0.003 
GeV, in which, ±0.016 is the error from theory (here, 
we mean that the upper and lower values of m-D^(o-) 
are determined by the upper and lower limits of s^-, 
as done in Ref.[ll|) and ±0.003 is the error from input 
parameters. We can see that the result of GSR is larger 
than the value of BSR[il|, 1.94 ± 0.03 GeV, and closer 
to the experimental value 1.969 GeV (s. Table HH). The 
values of /£)^(o-), '^'0^(0+) ^^d /d^(o+) are obtained in a 
similar way. The Numerical results are listed in Table HIl 
We can see from Table HIl cleanly that the results of GSRs 
are in accordance well with experiment than others. 

For checking the self-consistency of GSRs, we also com- 
pare the l.h.s. of Ea PH)) with the r.h.s., using the center 
values determined from GSRs. We have found that the 
two sides of GSRs, Eq. ([T0|. are compatible from each 
other very well. 



As a summary, in the framework of the Gaussian sum 
rules, the masses of the Ds{0^) mesons are investigated 
from a view-point of ordinary light-heavy system. The 
GSRs for the masses of Ds{0^) mesons and the cou- 
plings of Dg{0^) mesons to the corresponding currents 
are derived by means of the Laplacian transformation 
to the usual Borel sum rules. Using the standard input 
of QCD non-perturbative parameters, the corresponding 
mass spectra and couplings of the currents to the Ds{0^) 
mesons are obtained. By comparing both sides of the 
GSRs, we have shown that there exists a stability regions 
within which both sides of the sum rules are matched very 
well. Our results are mD^(o-) = 1-968 ± 0.016 ± 0.003 
GeV and mD,(o+) = 2.320 ± 0.014 ± 0.003 GeV, which 
are in accordance well with the experimental data, 1.969 
GeV and 2.317 GeV. Finally, it worth noting that the 
1^*^(2317) is treated as the lowest resonance of the 0+ 
cs meson in our calculation, and furthermore, its mass 
is even lower than that of the 0^ cd meson estimated in 
the similar way (this result will be published elsewhere). 
Therefore, -D*^(2317) may be considered to be the lowest 
0+ charmed meson in our treatment. 
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